Reaction dynamics of weakly-bound nuclei at near-barrier energies: impact of 
incomplete fusion on the angular distribution of direct alpha-production 
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The classical trajectory model with stochastic breakup for nuclear collision dynamics of weakly- 
bound nuclei is further developed. It allows a quantitative study of the importance of incomplete 
fusion dynamics in the angular distribution of direct alpha-production. Model calculations indicate 
that the incomplete fusion contribution diminishes with decreasing energy towards the Coulomb 
barrier, notably separating in angles from the contribution of no-capture breakup events. This 
should facilitate the experimental disentanglement of these competing reaction processes. 
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Introduction. Nuclear physics research has entered 
a new era with developments of radioactive nuclear 
beam facilities, where nuclear reactions are the pri- 
mary probe of the new physics, such as novel struc- 
tural changes. In those facilities, the low-energy nu- 
clear reactions research is highly focused on understand- 
ing astrophysically-important reaction rates involving ex- 
otic nuclei. These are often weakly-bound with a few- 
body, cluster structure that can easily be dissociated in 
their interaction with other nuclei. Understanding the 
breakup mechanism and its impact on nuclear reaction 
dynamics is essential. A major consequence of breakup 
is that a rich scenario of reaction pathways arises, such 
as events where (i) not all the resulting breakup frag- 
ments might be captured by the target, termed incom- 
plete fusion (icf), (ii) the entire projectile is captured by 
the target, called complete fusion (cf), and (iii) none of 
the breakup fragments are captured, termed no-capture 
breakup (ncbu). 

Since the availability of intense exotic beams is still 
limited, extensive experimental research has recently 
been carried out exploiting intense beams of stable 
weakly-bound nuclei, such as 6 ' 7 Li and 9 Be [H, 0, Un- 
derstanding the effect of their breakup on near-barrier 
fusion has been a key aspect of these investigations 
These have definitively demonstrated that breakup sup- 
presses the above-barrier fusion cross sections. Most 
recently, experimental activities are focused on disen- 
tangling breakup and competing reaction mechanisms 
from inclusive and exclusive coincidence measurements 
[E IE S B 9 • A challenge is to obtain a complete quan- 
titative understanding of the breakup mechanism and its 
relationship with near-barrier fusion. This research is 
gui ded by complete sub-barrier breakup measurements 

Theoretical works have addressed the low-energy re- 
action dynamics of weakly-bound nuclei using quantum 
mechanical, classical and mixed quantum-classical ap- 
proaches [ill, [H, EE Q EE EE H31 • Among these, the 
continuum-discretised coupled channels (cpcc) frame- 
work has been very successful [ill EE IxtI I - However, 
existing quantum models have limitations IE], as they 



cannot calculate integrated icf and CF cross sections un- 
ambiguously. Neither, after the formation of ICF prod- 
ucts, can they follow the evolution of surviving breakup 
fragment (s) since ICF results in depletion of the total few- 
body wave function. 

These difficulties are overcomed by the three- 
dimensional classical dynamical reaction model suggested 
in Ref. [2l| . It allows a quantitative description of ICF 
and CF of weakly-bound (two-body) projectiles, relating 
the sub-barrier ncbu to the above-barrier CF suppres- 
sion. 

I have further developed this approach. The key new 
aspect is the time propagation of the surviving breakup 
fragment and the ICF product, allowing the description 
of their asymptotic angular distribution and dynamical 
variables. This development should be very useful in 
(i) current experimental activities aimed at disentangling 
competing reaction mechanisms from asymptotic observ- 
ables such as alpha-production yields [9(, and (ii) appli- 
cations to 7 ray spectroscopy [IE HE] ■ The approach is 
quantitatively illustrated with a simplified test case. The 
methodology will be explained first. Afterwards, results 
are presented and discussed, and finally a summary is 
given. 

Methodology. Details of the approach are given in Ref. 
[2l| . and its implementation in the platypus code is 
described in Ref. [22]. Only the main features of the 
model are highlighted here: 

(i) The target T is initially at rest in the origin of the 

laboratory frame, whilst the weakly-bound (two- 
body) projectile P approaches the target (along the 
z-axis) with incident energy Eq and orbital angu- 
lar momentum L$. For each Lq (chosen to be an 
integer number of H) an ensemble of N incident pro- 
jectiles is considered. Including the P — T mutual 
Coulomb and nuclear forces, classical equations of 
motion determine an orbit with a definite distance 
of closest approach R min (E 07 L ). 

(ii) The complexity of the projectile dissociation is em- 
pirically encoded in a density of (local breakup) 

a function of the projectile- 



probability Vbu{R), 
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target separation R, such that V B ' u (R)dR is the 
probability of breakup in the interval R to R + dR 
(see Appendix A in [22]), A key feature is that for 
a given projectile-target combination, both experi- 
mental measurements |23j and CDCC calculations 
[2ll ] indicate that the integral of this breakup prob- 
ability density along a given classical orbit is an 
exponential function of its distance of closest ap- 
proach, R min (E ,L ): 

P B u(Rrmn) = 2 / V^ v (R)dR 

JRntin 

= A exp(-ai? mm ). (1) 

Consequently, V BU (R) has the same exponential 
form, V BU (R) oc exp(— aR). It is sampled to deter- 
mine the position of breakup in the orbit discussed 
in (i). In this position, the projectile is instanta- 
neously broken up into fragments Fl and F2. These 
interact with T, and with each other, through real 
central two-body potentials having Coulomb bar- 
riers V B 3 at separations Rg, i,j = 1, 2,T, i ^ j. 
The instantaneous dynamical variables of the ex- 
cited projectile at breakup, namely its total inter- 
nal energy £12, its angular momentum £12 and the 
separation of the fragments d\2 are all Monte Carlo 
sampled [2 11 ] . Having fixed the position and dynam- 
ical variables of the excited projectile fragments at 
the moment of breakup, the instantaneous veloc- 
ity of the particles Fl, F2 and T is determined by 
conservation of energy, linear momentum and angu- 
lar momentum in the overall center-of-mass frame 
(see Appendix B in [22j). These breakup initial 
conditions are transformed to the laboratory frame 
where the three bodies are propagated in time. The 
calculated trajectories of Fl, F2 and T determine 
the number of ICF, CF and ncbu events, fragment 
Fj being assumed to be captured if the classical tra- 
jectories take it within the fragment-target barrier 
radius R 3 B . 

(iii) From the N breakup events sampled for each projec- 
tile angular momentum Lq, the numbers of events 
Ni in which i = (ncbu), 1 (icf), or 2 (cf) frag- 
ments are captured determine the relative yields 
Pi = Ni/N of these three reaction processes af- 
ter breakup, with Pq + Pi + P2 = 1. The abso- 
lute probabilities Pi(E ,L ) of these processes are 
expressed in terms of the relative yields and the 
integrated breakup probability over the whole tra- 
jectory PBlj(Rmin): 

P()(Eo,Lo) = PBu(Rmin)Po, (2) 

Pi{E ,L ) = P BU (R mm )Pi, (3) 
P 2 (E 0l Lo) = [l-P B u(Rmin)}H{L cr -Lo) 

+ P B u(Rrmn)P2, (4) 



where H (x) is the Heaviside step function and L cr 
is the critical partial wave for projectile fusion. The 
cross sections are calculated using 

ai(Eo) =ttA 2 ^(2L + l)Pi(E , L ), (5) 

where A 2 = ti 2 /[2mpEo] and mp is the projectile 
mass. 

Beside the absolute cross sections ([5]), asymptotic ob- 
servables, such as the angle, kinetic energy and relative 
energy distributions of the fragments from NCBU events, 
are calculated by tracking their trajectories to a large 
distance from the target. 

Crucially, for the ICF events, the time propagation of 
the ICF product and the surviving breakup fragment is 
now incorporated into this picture. The captured frag- 
ment reaches the target radius forming the ICF product, 
whilst the other fragment flies away. At this moment, 
the three-body propagation turns into a two-body prop- 
agation, with definite interaction potentials and initial 
conditions. These are given by the position and veloc- 
ity of the three particles, at the moment when the ICF 
product is formed. 

This approach is here applied to the test reaction 8 Be 
_l_ zospj-, This is because a crucial model input is the 
local projectile breakup probability V BU (R), whose form 
deduced from PBu{Rmin) in U} is known thus far only 
for the reaction 9 Bc + 208 Pb from sub-barrier breakup 
measurements [23]. Very slightly changing the empirical 
V BU {R) from the 9 Be + 208 Pb experiment, the classical 
model CF and ICF excitation functions for 8 Be + 208 Pb 
very well agree with the experimental data for 9 Be + 
208 Pb at above-barrier energies (see Figs. 1 and 3 in [211]). 
Thus, the present outcomes for the angular distribu- 
tion of direct alpha-production are expected to represent 
those of the 9 Bc + 208 Pb reaction. The optimal breakup 
function PBu(Rmin) has parameters A — 5.98 x 10 3 and 
a = 0.85 fm -1 [see Eq. JTJ], The nuclear interaction 
between the alpha particle and the ICF product 212 Po is 
the Woods-Saxon (WS) potential well (V, r, a) = (33.98 
MeV, 1.48 fm, 0.63 fm) deduced from the global Broglia- 
Winther parametrization [24| . (Please note that in the 

1/3 

potential the radius parameter is multiplied by A T .) 
The rest of the model parameters are the same as in Ref. 

Results and discussion. Figure [1] shows the angular 
distribution of direct alpha-production for three incident, 
laboratory energies near the P — T s-wave Coulomb bar- 
rier (39.9 MeV), namely (a) E = 65, (b) 55, and (c) 
45 MeV. The contribution of the ICF and NCBU events is 
represented by thick solid and thick dashed lines, respec- 
tively. Their sum is represented by the thin solid line. 
It has a symmetric shape at energies well-above the bar- 
rier [panels (a) and (b)], adopting an asymmetric form 
as the incident energy decreases [panel (c)]. While the 
contribution of the ICF and ncbu events appears to be 
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TABLE I: The absolute ICF, ncbu and CF cross sections for 
the laboratory energies studied. 
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FIG. 1: Angular distribution of direct alpha-production for 
8 Be + 208 Pb for three laboratory energies Eo: (a) 65 MeV, 
(b) 55 MeV, and (c) 45 MeV. With decreasing energy towards 
the s-wave Coulomb barrier, the ncbu events dominate, sep- 
arating its centroid substantially from that of the ICF events. 
The total alpha-production distribution becomes asymmetric. 



similar at well-above barrier energies, the ncbu contribu- 
tion gradually dominates with decreasing energy towards 
the barrier. Here, its centroid significantly separates from 
the centroid of the ICF contribution. Both centroids shift 
to higher angles as the incident energy decreases, due 
to the reduction of relative partial waves affecting these 
reaction processes. However, the ncbu centroid always 
remains lower than the ICF centroid, as higher partial 
waves contribute to the ncbu process (see Fig. [5]). Fig. 
[5] presents the incident angular momentum distribution 
of the ICF (solid line) and NCBU (dashed line) processes 



E (MeV) 


aicF (mb) 


o ncbu (mb) 


ocf (mb) 


45 


84.7 


238.84 


149.2 


55 


252.3 


255.1 


624.6 


65 


273.6 


259.3 


981.5 



15 



10 



_Q 

E 








FIG. 2: ICF and ncbu cross sections as a function of the 
relative angular momenta Lo for 8 Be + 208 Pb at Eo = 45 
MeV. The contribution of high-partial waves shifts the NCBU 
distributions in Fig. [Tjto smaller angles, with respect to the 
ICF distributions. 



at a laboratory energy of Eq = 45 MeV. For complete- 
ness, the absolute icf, ncbu and CF cross sections are 
given in Table [I] for the laboratory energies studied. 

Summary. The classical trajectory model with 
stochastic breakup has been extended to the calculation 
of asymptotic observables associated with the ICF pro- 
cess. This development should be very useful for separat- 
ing the contribution of competing reaction mechanisms 
of weakly-bound nuclei from inclusive and exclusive mea- 
surements, such as the total alpha-production yield. This 
may also be affected by other direct processes, such as 
transfer [25| . which are not included in the model yet. 
The inclusion of transfer as well as the development of a 
unified quantum description remain as great theoretical 
challenges. The former may be carried out exploiting the 
concept of transfer function [26j], and the latter through 
a time-dependent density-matrix approach incorporating 
the concept of quantum decoherence [18|, [27| • Neverthe- 
less, for the first time, the impact of ICF dynamics on the 
angular distribution of direct alpha-production at near- 
barrier energies is quantitatively presented. The ICF con- 
tribution diminishes with decreasing energy towards the 
barrier, making the direct, total alpha-production dis- 
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tribution asymmetric. The ICF and ncbu contributions knowlcdged. 
are clearly separated in angles at energies close to the 
barrier, facilitating the experimental disentanglement of 
these competing reaction processes. 
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